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Abstract 

In this work the robustness of a simple cellular automaton devel- 
oped by Bassler and PaczuskiQ to describe the critical state in type- 
II superconductors is studied. Two different configurations of pinning 
centers are introduced and a new universality class is found. The nu- 
merical values of the critical exponents were calculated following two 
scaling techniques to ensure the validity of our results. 



1 Introduction 

In type-II superconductors, if the external magnetic field is greater than the 
first critical field H&, vortices penetrate the material. These vortices move 
freely in ideal superconductors if an electrical current is applied, destroying 
superconductivity ||. However, in real materials, the existence of pinning 
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centers prevent vortex motion if the current is less than a certain value, called 
critical current. 

Bean in 1966 |3j proposed that the distribution of magnetic field inside 
a real sample, when the external magnetic field increses over Hci, can be 
represented by a linear profile in a H(x) graph, where the slope corresponds 
to the critical current density. 

In 1989 Bak et al Q, in an attempt to explain the behavior of many dy- 
namical systems, developed the concept of self organized criticality (SOC). In 
their approach, after long time of avalanches evolution, those systems reach 
a complex state which exhibits a power law in their avalanche size distribu- 
tions and in their avalanche duration distributions, and f~@ noise. Sandpiles 
have become a paradigm of these sytems, and their closed similitude with the 
critical state of superconductors |J quickly called the attention of researchers 
in the field § §. 

In fact, Field et al in 1996 experimentally observed SOC in vortex 
avalanches at the inner wall of a low-T c superconducting hollow cylinder 
submitted to slowly ramped axial magnetic fields. In a similar experiment 
performed on a low T c thin film ring, Nowak et al flnj observed either poten- 
tial or peaked distributions of vortex avalanche size, depending on tempera- 
ture. Finally, in a recent experiment, Behnia et al [JTTj] measured the internal 
vortex avalanches in a low-T c film submitted to an increasing magnetic field, 
and found power distributions of avalanche size only for temperatures higher 
than 3K. 

On the other hand, many simulations had been performed on vortex 
avalanches. In 1994, Richardson et al fl2| developed a molecular dynamic 
model showing the discrete evolution of the magnetic profiles inside the super- 
conductors. Also using molecular dynamic techniques, Olson et al |13| tried 
to mimic Field's experiment. They found power distributions of avalanche 
sizes for high densities of pinning centers. However, none of these simula- 
tions was able to correctly estimate the critical exponents involved in the 
SOC theory for this system because of strongly finite size effects. 

More recently, Bassler and Paczuski |TJ proposed a simple cellular au- 
tomata to model the dynamics of magnetic flux motion at a length scale 
larger than the range of vortex interaction for a system with a random distri- 
bution of pinning centers. They obtained, using finite size scaling techniques, 
the existence of critical exponents for a wide distribution of parameters in 
their simulation, and hence, of SOC. In this work, we explore the robustness 
of this cellular automaton varying the distribution and depth of the pinning 
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centers. 

In the next section we explain Bassler and Paczuski's model. In section 3 
we present and discuss the results obtained using two novel pinning centers 
configurations: periodic distributions and combined periodic plus random 
distributions. Finally the conclusions are given. 

2 Model 

Bassler and Paczuski's cellular automaton [p]] is a two dimensional hexagonal 
lattice where each site is occupied by m(x) vortices. Vortices in the site x can 
move towards the site y only if the force acting on them, in that direction, 
is greater than zero. This force is calculated using the following formula: 

F = V(y) - V(x) + (m(x) - m(y) - 1) + 

r(m(xi) + m(x2) — m(t/i) — m(?/ 2 )) (1) 

where x\, x 2 and yx, y 2 are the other nearest neighbors of x and y, respec- 
tively, and V(x) and V(y) are the strengths of the pininng centers at those 
sites. If the distributions of pinning centers is random, V(x) = p, with a prob- 
ability q, and V(x) = with a probability 1 — q. The parameter r (where 
r < 1) characterizes the long distance action of the next nearest neighbours. 
If there is more than one unstable direction, one of them is chosen at random. 
All lattice site are updated in parallel, and at each site only one vortex can 
move on a particular update. Periodic boundary conditions are applied to 
the top and the bottom of the lattice. The vortices that reach the right edge 
of the system are removed, while they are not allowed to abandon the system 
through its left edge. 

An avalanche begins by randomly choosing a site at the left edge of the 
system and adding one vortex to it. It continues with the consecutive update 
of the lattice sites until no more unstable sites persist. Once the lattice 
is again stable, another vortex is added. The avalanche size is defined as 
the number of toplings corresponding to the addition of one vortex while 
the avalanche duration is defined as the number of updatings necessary to 
complete one avalanche. 

To characterize the system, in analogy with other SOC models, Bassler 
and Paczuski[|l| proposed and proved the following scaling ansatz for the 
probability distribution of avalanches sizes ([J), and avalanches durations @: 
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P(s,L) = s^g(s/L D ) 



(2) 



P(t,L)=t-^g(t/L z ) (3) 

and obtained the following set of scaling exponents for a wide range of pa- 
rameters values in the model r = 1.63 ±0.02, D = 2.7 ±0.1, r t = 2. 13 ±0.08 
and z = 1.5 ± 0.1 which are related by the scaling relations: 

t(2D - 1) = 1 (4) 
D(r-l) = z(r t -l) (5) 
In our work, to ensure the validity of our results, we also used the fol- 



lowing scaling ansatzs |TJ] for the avalanches sizes and avalanches durations, 
respectively: 

P(s,L) = L-?f(s/L D ) (6) 

P(t,L) = L-f(t/L z ) (7) 

Therefore two new scaling exponents characterize our system j3 and to. 
They are related to the previous ones by: (3 = tD and uj = r t z, and satisfy 
the following scaling relations: 

P = 2D - 1 (8) 

and 

UJ = f3 + z-D (9) 

With this notation, Bassler and Paczuski's results U give (3 = 4.4, D = 2.7, 
u = 3.19 and z = 1.5. 



3 Results and Discussion 

Following Bassler and Paczuski|l|], we first assumed a random distribution of 
pinning centers, with parameters (r,p,q) = (0.1; 3; 0.1). In figure 1 appears 
the collapse of our results for dimensions L = 100, 160, 200 using the scaling 
ansatz (|6|) and ([?]). The critical exponents for these conditions are in good 
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agreement with reference [|I[] . Similar results were obtained for other param- 
eters within the range already checked in 0, in all case after a 10 7 m.c.s 

Then, we tested the model using a combination of periodic and random 
pinning distributions. This kind of configuration introduces a new parameter 
to geometrically characterize the system i.e, the distance between the peri- 
odic pinning centers a. In our simulations we used r = 0.1, and a constant 
density of random pinning centers q = 0.1. Like real irradiated superconduct- 
ing materials, we used different pinning strenghts p\ and P2 for periodic and 
random pinning centers, respectively. The periodic pinning centers were al- 
ways chosen stronger or equal to the random ones trying to mimic irradiation 
effects. 

For the sets of parameters (puP2',o,) = (10;1;20), (10; 1; 10), (5; 1; 10), 
(10; 1; 4) and (20; 1; 4) the critical exponents conserved the previously re- 
ported values of: t = 1.63±0.04 , Tt = 2.13±0.09, (3 = 4.4±0.1, D = 2.7±0.1, 
uj = 3.2 ± 0.2 and z — 1.5 ± 0.1 for different system dimensions, indicating 
the robustness of the system, (see the two data collapse presented in figure 
2). 

We also performed simulations using only a periodic pinning configura- 
tion. It was demostrated that the system also conserved the critical expo- 
nents for the sets of parameters (p; a) = (5; 10) (5; 20), (5; 4), (1; 10) as is 
shown in figure 3. However, as the data collapse in figure 4 displays, us- 
ing lower values of a and stronger pinning centers (for example: (10,4) and 
(20, 4)) caused a decreased in the exponents j3, r, r t and to to (3 — 3.2 ± 0.1, 
r = 1.45 ± 0.02, r t = 1.7 ± 0.08 and u = 2.6 ± 0.2, while z = 1.6 ± 0.1 
and D = 2.2 ± 0.1, in good agreement with the scaling relations (f|) (|5]) ([8]) 
and @.This suggests the existence of a new universality class for dense and 
strong periodic pinning configurations. 

To study the origin of this new universality class we made calculation 
using a superposition of random and periodic pinning configurations with 
equal potentials p± = p2- The scalings results, presented in the figure 5, show 
that systems with parameters (pi',P2',a) = (10; 10; 4) conserved Bassler and 
Paczuski's|l]] exponents values. 

These results and those obtained using only a periodic configuration with 
the set of parameters (p; a) = (5; 4) (see figure 3) proved that the new uni- 
versality class appears due to the joint action of the strong periodic pinning 
centers and the neglegible influence of the random pinning distribution. If in 
the system one of these conditions is absent, it will have a "normal" evolution. 
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The existence of the new universality class could be explained based on 
the following qualitative argument. The presence of strong and correlated 
pinning, the neglegible influence of the random pinning and the low values 



of a could produce "magnetic traps" [|13|] were vortices oscillate between 
neigbour pinning centers, increasing the number of bigger avalanches, and 
then reducing the characteristic exponents, r, r t , (3 and lu. 

Finally an experiment is suggested to check the existence of this new uni- 
versality class. It consist in performing an avalanche detection experiment 



similar to the one reported in reference [11|, but before an after irradiating 



the sample with heavy ions in order to obtain arrays of periodic " strong" pin- 



ning centers as reported by Harada et al [15|. If some relevant experimental 



parameters such as the magnetic field sweep rate and the temperature are 
conveniently tuned, critical exponents close to those reported in [l] should 
be achived before irradiation. After the irradiation (under the same set of 
parameters) the combination of "weak" random pinning plus "strong" peri- 
odic pinning is expected to move the critical exponents towards the values 
reported here. 



4 Conclusions 

A cellular automaton model for vortex avalanches was analized in the pres- 
ence of periodic-random and periodic distributions of pinning centers. The 
robustness of the model was proved for a wide range of parameters and the 
existence of a new universality class for strong and dense periodic distribu- 
tion of pinning centers was shown. A qualitative argument to support the 
reasons for the appareance of this new universality class were given. 
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Figure captions 



Figure 1 Finite scaling plot using (^) L = 100,160,200 and (r;p;q) = 
(0.1; 3; 0.1). Inset: finite scaling plot using (0) 

Figure 2 Finite scaling plot using (|7|) L = 100,160,200 and (r;pi;p2', a) = 
(0.1; 10; 1; 10). Inset: finite scaling plot using (§) 

Figure 3 Finite scaling plot using (|6|) L — 100,160,200 and (p; a) = (5; 4). 
Inset: finite scaling plot using (|3|) 

Figure 4 Finite scaling plot using (0) L = 100, 160, 200 and (p; a) = (10; 4). 
Inset: finite scaling plot using (||) of 

Figure 5 Finite scaling plot using (|7|) L = 100,160,200 and (r;p 1 ;p 2 ;a) = 
(0.1; 10; 10; 4). Inset: finite scaling plot using (Q) 
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